16.451 Lecture 17: Beta Decay Spectrum ( Ref:F&Hchll ) Nov. 3, 2005 


n — » p + e + v e 

(and related processes...) 


Goals: 


• understand the shape of the energy spectrum 

• total decay rate sheds light on the underlying 

weak interaction mechanism 



Starting point: "Fermi's Golden Rule" again! (lecture 6) 


matrix element 



(transitions / sec) 


M if = jVf V(r) y/i d 3 r 


Pf - dnt dEr 

density of states 



Matrix element M if ? 


Simplest model is to take a pointlike interaction with an overall energy scale " G 

(Fermi, 1934 - almost right!) 



My = G J y* p f (r ) i/>* (r ) tj>* v (r ) i// n i (r ) d 3 r 


(the interaction is proportional to the wavefunction 
overlap of initial and final state particles at the same 
point in space) 



Standard Model description: an extended 
interaction, but the range is only about 0.002 fm 
which is just about zero! 

The Standard Model can 'predict' the value of G 
in terms of model parameters, whereas in Fermi's 
theory, it remains to be determined from experiment. 



Spin considerations: electron and neutrino 3 

There are two possibilities for the angular momentum coupling of the two leptons! 

S e + S v = V => s = 0 or 1 

For neutron decay: n —> p + e + V e 

Angular momentum: — = A + S can contribute to neutron decay! 


Subtle point: because the leptons are emitted with a def inite helicity 

we can deduce a correlation between their directions of motion in the two cases: 


Fermi Decay, S = 0: 

s e (. LH ) 

► * — 

◄ ◄ 

Sy (J ?») 


e- and v travel in the same direction 


Samow-Teller decay. 

S = 1: 

s e (LH) 

Pe 



Sy (. RH ) 

Pv 

e- and v travel in opposite directions 



More on spin correlations: 


Fermi Decay, S = 0: n —> p + e + V e 

Leptons travel in the same direction:: 


p p 


Recoiling proton spin is in the same direction 
as the initial neutron spin. 


Gamow - Teller Decay, S = 1 : 

Leptons travel in the opposite direction: 


p p 


Recoiling proton spin is in the opposite direction 
as the initial neutron spin, i.e. a "spin flip" 


S = 1 







How to proceed? 5 


As before, assume a pointlike interaction, but allow for different coupling constants 
for the Fermi (F) and Gamow-Teller (GT) cases. 

Fermi case, S = 0: (coupling constant: "G v " because the operator transforms 

like a space vector.) 

M if = Mp = Gy JV* >f (r)0*(r)0*(r)tr n 'i(r) d 3 r 

Gamow-Teller, S = 1: (coupling constant: n Q A " because the operator transforms 

like an axial vector, i.e. like angular momentum.) 

M, f = M gt = G A jy/* p f (r)0*(r)f v (r)y/ ni (r) d 3 r 

Experimentally, the coupling constants are very similar: 

G a /G v = - 1.25 

(These are evaluated by comparing different nuclear beta decay transitions, where 
angular momentum conservation restricts the total lepton spin states that contribute) 



Transition rate: 


n ^ p + e + v e 


(can proceed with S = 0 or 1) 




for the neutron: 


r— d • 

-"■h 

II 

^ £ 

M if 

P f ~ (Gy +3 G\) 



since there are 3 times as many ways for the leptons to be emitted with 5 = 1 
(m s = 1, 0, -1) as with 5 = 0. 

For now, let us work out a generic matrix element, since the expressions are the 
same for both apart from the coupling constants: 

M if = G d 3 r 


Pr 

i 

recoil 



Electron and neutrino are represented by 
plane wave functions of def inite momentum: 



ip*r /h 

Vv 


etc. 



Matrix element: 


7 



o N 


recoil v n, 

_ S A 

<7 neutrino 


M if = G jv/* p j(r)^*(r)£(r)i// n j (r) d 3 r 


Pr = ~(p + q ) 






[ p R .r/h 


The integral for M if extends over all space regions for which the nucleon wave functions 
(n,p) are non-zero: R max ~ 1 fm (in nuclei, use R ~ 1.2 A 1/3 fm ) ... 

But, the recoil momentum p R is no larger than the Q-value for the reaction, ~ MeV ... 


p R *r/h < 


Q R tmx 

Tic 


1 MeV.fm 

« 1 

197 MeV.fm 



This is a great simplification: the lepton wave functions are just a constant over the 
region of space that matters to calculate the matrix element © 



Summary so far for the transition rate calculation: 


o 2^- 


2 

271 

G If* 7 3 

2 

» 

M{f 

Pf 

fl 

y2 |J Vf.p Yi,n d r 

Pf 


The remaining integral is referred to as the nuclear matrix element: 



When beta decay occurs in a nucleus, the initial and final wave functions of the proton 
and neutron need not be exactly the same, so in general: 



nuclear 


< 1 


However, in the case of the free neutron, there are no complicated nuclear structure 
effects, and so the matrix element is identically 1: 




2x_ C?_ 

h y 2 



When this occurs in a nucleus, the beta decay 
rate is the fastest possible, and the 
transition is classif ied as "superallowed" 



Density of states factor: 


Just like the calculation we did for electron scattering, but now there are two light 
particles in the final state! 



o. 


recoil v x 

Q neutrino 


We want to work out 
the number of equivalent 
final states within energy 
interval dE f of E f . 


Final state momenta are quantized in volume V (lecture 6) 


dn = dn e x dn v 


r 2 V ^ 

An p dp 

V h*J 


x 


( 2 

An q dq^r 

V h 6 J 


Ef = Me 2 + me 2 + K r + K e + cq 
But the nucleon is much heavier than the other particles: K R = ( p R ) 




K e =const 



9 


dn 

dE f 


I2M = 0 




Put this all together for the density of states factor: 


10 



And finally, for the transition rate: 



Exercise: plug in all the units and check that the transition rate is in sec 1 


Notice: This is actually a partial decay rate, because the electron momentum p is 
specif ied explicitly. X if here gives the rate at which the decay occurs for a given 
electron momentum falling within dp of p -» this predicts the momentum spectrum! 



Electron momentum spectrum: 


N(p ) dp = N 0 A if 
=> N(p) 


(i const.)x p 2 q 2 dp 
(const.) x p 2 q 2 = 


approx: K R 

. / 


= 0 


(const) x p 2 (Q-K e ) 2 


Pr 

I 

recoil 



Predicted spectral shapes, 
Krone, f igure 9.2: 

(plotted for Q = 2.5 MeV, 
not the neutron!) 

(Note: max. K(e ) = Q ) 




Comparison to Reality: (Krane, fig. 9.3: e+ and e- decays of 64 Cu) 12 


Too many low energy e - 



Too few low energy e + 



Discrepancy: neglect of Coulomb effects in the final state. 


Key point: Coulomb distortions of the energy spectra occur AFTER the electron/positron 
are emitted in the weak decay process. 


Modified density of electron/positron states: 



"Fermi function", depends on 
the charge 71 of the "daughter 
nucleus" (final state) and the 
electron/positron momentum 


Approximate correction factor for p± decay: 


F ± (Z\p ) = 


x 


1- e 


- x 


__ 2n a Z' v e l 1 

x = + , p = —, a = = 

P c 47rs 0 hc 137 



Modified electron/positron spectrum prediction: 


N(p) 


C p 2 ( Q-K e f F ± (Z',p), 


C = 






Fermi-Kurie Plot 


14 


Idea: for "allowed decays", corresponding to our approximation: e ip R -r/h = x 

inside the nucleus, the electron energy spectrum can be "linearized" if one 
accounts for the Coulomb distortion via the Fermi function F(Z',p): 



66 Ga -> 66 Zn decay, Phys. Rev. 129, 1782 



~ Q ~ K e 

-V 

linear function, 
endpoint Q 


K e + m e 





Neutrino Mass effect: 3 H — » 3 He + e + v e (Q = 18.6 keV) 


15 


Idea: shape of the electron energy spectrum near the endpoint (Q) is sensitive to the 
mass of the electron antineutrino: 


Pr 

I o 

recoil 




neutrino 


recall: Q = K R + K e + K y 


WhenK e = Q, K R = K v -> 0. if m v ^ 0, then 
in this limit, mass effects are most pronounced. 



Figure 2: The electron energy spectrum of tritium (3 decay: (a) complete and (b) narrow region 
around endpoint Eq. 1 he 13 spectrum is shown for neutrino musses of 0 and I eV. 


figure from "KATRIN" proposal, 2001: (Karlsruhe Tritium Neutrino expt.) 
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Figure 4: Results of tritium $ fleetly experiments on the observable fflQ over the IelsL decade. 

• Best direct upper limit: m v < 2.2 eV 

• from Sudbury neutrino observatory and other experiments, we have convincing 

indirect evidence of nonzero neutrino mass that is much smaller than this 



Web sites: http://cupp.oulu.fi/neutrino/nd-mass.html . http://www-ik.fzk.de/~katrin/ 



